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A physical model is proposed for the estimation of the screech amplitude in 
underexpanded supersonic jets. The model is based on the hypothesis that the 
interaction of a plane acoustic wave with stationary shock waves provides 
amplification of the transmitted acoustic wave upon traversing the shock. 
Powell’s discrete source model for screech incorporating a stationary array of 
acoustic monopoles is extended to accommodate variable source strength. The 
proposed model reveals that the acoustic sources are of increasing strength with 
downstream distance. It is shown that the screech amplitude increases with the 
fully expanded jet Mach number. Comparisons of predicted screech amplitude 
with available test data show satisfactory agreement. The effect of variable source 
strength on the directivity of the fundamental (first harmonic, lowest frequency 
mode) and the second harmonic (overtone) is found to be unimportant with regard 
to the principal lobe (main or major lobe) of considerable relative strength, and is 
appreciable only in the secondary or minor lobes (of relatively weaker strength). 
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I. Introduction 


Supersonic jet screech constitutes an important technological challenge in the theory of jet 
noise. Screech refers to high-amplitude discrete tones superposed on the otherwise broadband 
acoustic spectrum observed in imperfectly expanded supersonic jets (Fig. 1). The shock noise 
produced by imperfectly expanded (choked) jets is very complex, and consists of two separate 
components, namely broadband shock noise (arising from shock-turbulence interaction) and 
screech tones. This shock noise is in addition to the broadband turbulent mixing noise in the jet 
[1-3]. In the nearfield, the very high intensity of these screech tones (as high as 170 to 180 dB) 
with a highly directional character with significant upstream propagating component can induce 
fatigue and cause serious structural damage and failure of aerodynamic vehicles in flight. A 
fundamental understanding of jet screech amplitude is thus requisite to the design of a new 
generation of rocket engines [4-6], turbo-jet aircraft engines, and to a number of other flow 
circumstances such as jet impingement and cavity resonance involving feedback cycles. 

Owing to pressure mismatch (nozzle exit pressure exceeds the ambient pressure), 
underexpanded supersonic jet structures manifest themselves in shock-cell structures involving 
normal/oblique shocks and regions of expansion (Fig. 2). Under certain conditions, a Mach disc 
is also formed. These shock-expansion units interact with instability waves, vortices, turbulence, 
and other stream disturbances in the viscous shear layer that surrounds the inviscid region. 

Powell [7-9] first identified the mechanism of supersonic jet screech in terms of a resonant 
feedback loop or cycle that is self-sustaining and consists of the following key processes (Fig. 2): 
downstream passage of flow disturbances in the jet shear layer initiated near the nozzle lip region 
(where the relatively thin shear layer is susceptible to instability); interaction of these stream 
disturbances with the nearly periodic and stationary oblique shock-cell structure of the jet 
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resulting in their amplification, with the ultimate generation and emission of intense acoustic 
waves or energy; upstream propagation of the acoustic waves immediately outside the jet shear 
layer (acoustic feedback); excitation and maintenance of stream disturbances in the vicinity of 
the nozzle lip by the upstream propagating acoustic waves (so-called receptivity), which then 
propagate downstream, thus closing the feedback cycle without any external interference. 
Experimental data, such as flow visualization measurements by Kaji & Nishijima [10], suggest 
that the shock-noise is most intense somewhere between the third and fourth shock-cell 
structures. For a detailed discussion on screech, Seiner [11] may be consulted. 

Powell’s theory [7] of stationary array of discrete monopole sources distributed through the 
shock-cells provided directivity patterns of screech that agree with measurements, as reported by 
Norum [12] for unheated jets. It also furnished a framework for the estimation of screech 
frequency based on shock-cell spacing. Fisher & Morfey [13] and Tam et al. [14] provided 
expressions for the prediction of screech frequency over a wide range of fully expanded jet Mach 
number Mj and jet total temperature T 0 . Tam’s expression for the screech is predicated on the 

so-called weakest-link theory, according to which the feedback acoustic waves represent the 
weakest link and control the screech frequency. Massey & Ahuja [15] expanded on Tam’s 
formulation and provided expressions for screech frequency for hot jets with different screech 
modes. The screech mode changes from axisymmetric to helical and flapping modes as the Mach 
number Mj increases. 

On the basis of Powell’s model of stationary point sources, Harper-Bourne and Fisher [16] 
experimentally and theoretically treated broadband shock noise (radiated to the farfield) by 
considering shock-turbulence interaction. Harper-Bourne and Fisher considered spatial 
coherence of the disturbance field to account for the broadband sound radiation to the farfield. 
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Generally speaking, the intensity of broadband shock-associated noise is more prevalent in the 
forward arc, and is a function only of the pressure ratio. A stochastic model for broadband shock 
noise was reported by Tam [17] who considered the nonlinear interaction between jet turbulence 
and stationary shock-cell system. Shock-turbulence interaction (shear layer and an oblique 
compression wave) was numerically investigated by Lui & Lele [18] including the dynamic 
motion of the shock during interaction with shear layer vortices. 

With regard to the screech process, our current understanding of the highly nonlinear flow 
dynamics and of the rate of amplification of stream disturbances is rather incomplete. Whereas 
screech frequency is controlled by the feedback acoustic waves, the screech intensity is governed 
primarily by the characteristics of the downstream propagating flow disturbances and their 
interaction with the shock-cell structures. Several detailed measurements [10, 12, 15, 19-23] and 
numerical viscous simulations [10, 24-29] have been reported and served to throw light on the 
subject. The experimental description of A-D mode switch (A u A 2 ,B,C,D) according to the jet 
operating parameters represents a notable step [12]. Despite the theoretical advances in screech 
frequency predictions, directivity character, modulations of instability waves by shocks [30-31] 
and unsteady shock motions and their role in noise generation, screech amplitude predictions 
remain elusive [32-33]. Owing to strong nonlinearity of the feedback loop, even an empirical 
formula for screech intensity is unavailable [24], 

The preponderating questions concerning the interaction of shock waves with stream 
disturbances and the excitation of disturbances by upstream-traveling acoustic waves largely 
remain unanswered. The purpose of the present work is to investigate the amplitude of screech 
based on a physical model for the interaction of stream disturbances (primarily regarded as 
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finite-amplitude acoustic waves) with shock waves. This paper is primarily based on the author’s 
recent work [6]. 

II. Proposed Physical Model 

A. Source Distribution 

The present model builds on Powell’s conception of stationary discrete acoustic monopoles 
(Fig. 3). However, in the present context account is taken of the variable strength of the 
monopoles, as the flow disturbances are amplified on their passage through the shock-cell 
system. As with Harper-Bourne and Fisher [16] the acoustic sources are regarded to be 
distributed along the lip line at locations in the jet shear layer where shock reflections take place. 
It is assumed that the sources S l ,S 2 ,S i are equally spaced. It is also assumed that relative to the 
source S 2 , the source S, leads in phase r and S 2 lags in an equal amount of phase r . 

The downstream-propagating stream disturbances near the supersonic side of the shear layer 
and the upstream-traveling acoustic disturbances on the subsonic outer edge of the shear layer 
form two key components of the feedback process. Attention is focused here on the amplification 
of downstream-traveling stream disturbances. 

B. Shock-Acoustic Wave Interaction 

1. Transmission of an Acoustic Wave through a Normal Shock 

The proposed model for the screech amplitude is based on the doctrine that the downstream- 
propagating stream disturbances in the screech cycle are. regarded as acoustic-like (nearly 
isentropic pressure fluctuations) and are of finite amplitude. The amplification of these 
disturbances during their passage is a consequence of the interaction between an acoustic wave 
(of finite amplitude) and a stationary shock wave. If the screech cycle is regarded as analogous to 
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a thermodynamic cycle (with energy transfers occasioned during shock interaction and scattering 
at the nozzle lip), the acoustic waves form the working medium for the screech cycle. Thus, it is 
likely that the upstream-propagating acoustic waves and downstream-traveling acoustic waves 
constitute the screech cycle. The scattering of upstream-propagating acoustic waves by the 
nozzle lip ensures that this is so, even though some of the incoming acoustic energy is 
transmitted to the hydrodynamic disturbances or instability waves in the shear layer in the 
receptivity process. 

Unsteady phase-averaged measurements of nearfield pressure fluctuations in a round jet at 
Mj =1.19 by Panda [21] indeed reveal that the downstream-propagating disturbances contain 

both hydrodynamic pulsations and acoustic waves, whereas upstream-propagating components 
are only acoustic [33]. 

When an acoustic wave (traversing with sound speed relative to the moving fluid) is incident 
upon the shock from ahead of it (Fig. 4), no reflected acoustic wave can arise because the flow 
ahead is supersonic, and the acoustic wave (longitudinal isentropic sound wave) is transmitted 
(refracted) downstream of the shock. At the same time, perturbations in entropy and vorticity 
(both are nonisentropic and nonacoustic; these isobaric, transverse waves are termed entropy- 
vorticity waves after Carrier) are generated, which are transported along with the flow behind the 
shock [34-35]. On the other hand, when an acoustic wave is incident upon the shock from behind 
it, no disturbances can arise in the medium ahead of the shock. The acoustic waves are merely 
reflected, and the entropy and vorticity waves are generated behind the shock. Thus, the 
refraction of the incident acoustic wave is assured as long as the approach flow is supersonic. 
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The main relationships developed here are based on transmission through normal shocks, 
although the shocks in the jet are not wholly normal but oblique. The Mach number normal to 
the shock M\ n is related to the approach Mach number M\ by the relation 

M\ n = M\ sin P\ 

where represents the shock wave angle (angle between the approach velocity and the plane of 
the shock). However, it is believed that the relationships based on the normal-shock assumption 
provide an upper limit to the screech tone intensity. 

The passage and linear interaction of plane disturbances (sound waves) with a shock wave 
upon normal incidence for an ideal gas was first treated by Blokhitsev [36] and Burgers [37]. 
When the perturbing waves are normally incident upon the shock, the transmission coefficient is 
expressed by [36-37] 


§h_ _ + 2 (r + IX + (3 r -W + I -r 

4>, (r + iXi + < + 2^,0 


where Sp { and Sp 2 denote the strength of the incident and refracted acoustic waves, 
respectively, M [n the normal Mach number of the flow upstream of the shock, M ln the normal 
Mach number of the flow downstream of the shock, and y the specific heat ratio. The 
downstream Mach number is obtained from the Rankine-Hugoniot normal shock relation [38] 


U 2 M| 2 „+2/Q-i) 

2 " [2y/(r-l)K,-l 
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Fig. 5 presents the variation of the pressure amplification fijpj ! &P\ as a function of incoming 
Mach number Mj. The pressure amplification is seen to rapidly grow as the upstream Mach 
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number is increased. For large values of Mi , the pressure ratio varies asMi . For oblique shocks 
(P < 90 deg), the pressure amplification is less than that for the normal shock. 

The result of the linear shock-acoustic wave interaction theory has been shown to be in close 
agreement with the numerical nonlinear Euler simulations by Zang et a. [39] for incident angles 
within about 30 deg of the critical angle (zero degree incident angle refers to normal incidence) 
for the acoustic as well as vorticity wave. Although linear theory is employed here for shock- 
acoustic wave interaction, the successive amplification is considered nonlinear. 

2. Successive Amplification (Shock-Cell Interaction) 

Based on the general behavior of the Mach number decay for typical underexpanded jet 
conditions, the following distribution of the approach (or upstream) Mach number for each tip of 
the shock-expansion unit (shock-cell reflection occurring in the shear layer) may be considered 
as a rough approximation: 

M\i = M r -{M r -l/sin/?i)(/-l)/(m-l) i = \,2,...m (2a) 

where M r is the upstream Mach number for the first shock-cell, M (/ - the incident Mach number 
in the / th shock-cell, and m the total number of active shock-cells under consideration. Thus, 
Eq. (2a) ensures that M (/ « M r for the first shock-cell (i = 1 ), and M Vl = 1/sin for the last 

cell ( i = m ). Based on the available data, the quantity M r is approximated as 

M r w 1 AM j (2b) 
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The quantity My could be more accurately obtained by a direct calculation, using numerical 
methods. In the present work, the choice of Eq. (2a) for approximating My is based on the 
consideration that it facilitates a simpler formulation of the model. 

The fully expanded Mach number M : is related to the pressure ratio p t / p a and the Prandtl- 

Glauert parameter (5 by 


M j = 

7 7 - 1 


( Pt 


VP a 


-1 


(2c) 


and 




(2d) 


with Pt and p a denoting jet total pressure and the ambient pressure, respectively. The Prandtl- 
Glauert parameter /? is thus related to the strength of the normal shock wave. The number of 
shock reflections for a convergent-divergent nozzle strongly depends on the nozzle design Mach 
number Mj and the operating pressure ratio. 

The relative source strength in each shock-cell can be determined from Eqs. (1) and (2). Eqs. 
(la) and (2a) assert that there is successive amplification of the sound waves during their passage 
through the shock-cell system, with the highest amplification rate occasioned in the first shock- 
cell. The overall (or total) amplification becomes 


4 Pm _ 



( o A 

vPm 

APin 

^$Pin A ) 


VpPm - 1 J 


( 3 ) 


where Ap in refers to the disturbance level upstream of the first shock-cell in the vicinity of the 
nozzle exit. 
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III. Screech Character 


A. Screech Amplitude 

The screech amplitude (so-called excess tone level) relative to the broadband spectrum SPL s 
is expressed as 

SPL s = 20 log ! o (A p m l A p in ) (4) 

where tsp m / Ap jn is determined from Eq. (3). Thus, Eqs. (la) to (4) enable us to compute the 

excess (or relative) screech amplitude as a function of jet Mach number. The excess tone noise 
represents the difference between the screech amplitude and the broadband level (both shock- 
induced and turbulence-induced) at the screech frequency in a given spectrum [15]. 

The above theory does not display a separate effect of jet temperature on screech tone 
intensity. Data of Massey et al [15] suggest that the tone intensity reduces with an increase in 
temperature. According to the measurements of Massey and Ahuja [40] mode switching occurs 
at elevated jet temperature, thereby affecting screech intensity. 

B. Screech Directivity 

In the present formulation, an extension of Powell’s treatment [7-9] is considered for the 
directivity of screech based on stationary point sources along a line; see Lord Rayleigh [41]. The 
farfield radiation is determined synthetically by combining these sources. With each source of 
strength of An and 4 7tH for the fundamental and the first harmonic respectively, the 
distribution of velocity potential for a point source is of the form 

$ = l e ik ( ct ~ r ) + !L e ik{ct-r) ( 5 ) 

r r 
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where c denotes the sound speed, r the radial distance, t the time elapsed, and k the wave 
number, which is related to the wavelength X of the feedback acoustic wave by 

X = 2 nlk (6) 

Following Powell, we assume that the phases rj = r 2 = r . For small ratios ois/r (where 
s denotes the shock-cell spacing), the potential for the fundamental of a three-cell shock system 
(for example) can be expressed as 

t / J l + ° 2 +° 3 ) D/Cosk (c l - r -S,) ( 7 ) 

The quantities aj and <33 respectively refer to the relative strengths of the second and the third 
shock-cell with the strength of the first shock-cell taken as unity. The phasing S t between 
consecutive monopoles is determined by shock-cell spacing and disturbance convection Mach 
number M c [ 12 ]. 

The directivity D f can be expressed as 

Df= 7 r \p2 + (l + <33 )cos k{cr - s cos #)] (8) 

(1 + + °3 ) 

with t = s / u c = s /( cM c ) ( 9 ) 

where 6 denotes the angle from the downstream flow direction and u c and M c refer to 

convective velocity and convective Mach number of the disturbance, respectively. Eq. (8) can 
be recast as 
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f (l + cr 2 +a 3) 


a 2 + (l + a 3 )cos — cos 0) 


XM„ 


Thus, the directionality of screech strongly depends on the shock cell spacing and the 
disturbance Mach number [7]. The quantity (l-A/ c cos#) is a Doppler factor incorporating the 
variation in retarded time and source phasing [16]. An expression for the directivity of the 
harmonic D h is provided by Eq. (10) with in replaced by4/r . In the special case ofa 2 = o 3 = 1 , 

Eq. (10) reduces to that derived by Powell [7] for uniform source strength. Eq. (10) can be 
extended to four or more sources. 

The shock-cell spacing s is obtained from experimental data and is related to the pressure 
ratio parameter /? defined by 


s/dj =f{p) 


For example, data of Harper-Bourne and Fisher [16] suggest that 


s/dj « 1.25/? 


Experimental data [7] suggest that the shock-cell spacing is related to the feedback acoustic wave 
length as 


s/X»0A (11) 

where X is an integer. 

Detailed measurements by Harper-Bourne [42] suggest that the shock-cell spacing s is not 
uniform, but varies somewhat along the jet. A non-uniform shock-cell spacing, which may affect 
the directivity, is unaccounted for in the present theory. Experimental investigations of Massey et 
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al. [15] also indicate that the effect of heating on the directivity of screech is largely unimportant. 
No direct effect of jet temperature on the directivity of screech is however evident from Eq. (10). 

C. Existing Relations for Screech Frequency 

Tam et al. [14] proposed for the screech frequency a semi-empirical expression of the form 

fdj/Uj=f{Mj, r ,T 0 /T a ) (12a) 

where T a represents the ambient temperature. The fully expanded jet diameter dj is related to 
the jet exit diameter and the nozzle design Mach number M d as 

dj/d = f(Mj,Md,r) ( 12 b) 

This formula agrees well with the data and accommodates the effect of jet temperature (hot and 
cold jets). Massey and Ahuja [40] provided separate expressions, similar in form to Eq. (12), for 
different screech modes depending on the Mach number range. 

The screech frequency increases with temperature since at a fixed Mach number, the 
convective velocity u c of the instability wave increases with temperature [15], The frequency is 

inversely proportional to the total time of travel in the screech cycle, (7j + T 2 ) , where 7] refers 
to the time for the downstream-propagating instability wave to traverse the shock-cell system 
starting from the nozzle exit, and T*i the time elapsed for the upstream-propagating feedback 
acoustic wave to reach the nozzle lip region. Thus, as the temperature increases, the value of 
T\ increases so that the frequency decreases. 
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IV. Results and Discussion 


In all the results presented here, except for the screech amplitude, all the sound pressure levels 
(dB) are taken relative to 20 n Pa. The quoted levels are in dB linear (unweighted). 

A. Streamwise Distribution of Acoustic Source Strength 

The streamwise distribution of the acoustic source pressure, as predicted by the present 
model for the normal shock case (y9j = 90 deg), is depicted in Fig. 6a for various Mach numbers. 
The successive amplification of the stream disturbances along the shock cells is evident, with the 
source strength increasing with the Mach number. The corresponding results for an oblique 
shock with = 60 deg are presented in Fig. 6b. The amplification rates for the oblique shock 
case are less than those for the normal shock case, as evident from Fig. 5. 

A comparison of the predicted amplification rates (for (5\ = 90 deg and = 60 deg) with the 
data of Kaji and Nishijima [10] for Mj = 1.64 and of Panda et al. [19] for Mj = 1.19and 1.42 is 

illustrated in Fig. 6c. The data are deduced from the original measurements in order to provide a 
direct comparison. In this context, it would be instructive to briefly present the data in the 
original form. 

The streamwise amplitude measurements by Kaji and Nishijima [10] for =1.64 (Fig. 7a) 

suggest that the highest intensity is located in the shear layer of the third or the fourth shock cell, 
whereas inside the jet the acoustic intensity is low. Interestingly, the data reveal low-pressure 
dips (dimple spots) in the source strength. According to Kaji and Nishijima [10] these dips are 
believed to be generated by the interaction between sources. There is a possibility that they could 
also be due to the expansion regions in the flow. 
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Detailed measurements by Panda et al. [19] for a circular jet over a wide range of Mach 
number ( M } = 1.19 to 1.51) show similar trends for the axial distribution of source strength, 

with the peak values as a whole showing an increasing trend with Mj , ultimately approaching a 

level of saturation followed by a relatively rapid decline in the screech levels (Fig. 7b). 

The foregoing data for the amplification rate are deduced here such that the relative source 
strength from the data at the first shock cell is taken as that given by the theory. Fig. 6c suggests 
that the proposed theory based on the normal shock for the acoustic source strength 
at Mj =1.2 is in close agreement with the data of Panda et al. [19] for all the shock cells. 

However, at Mach numbers of 1.4 and 1.6, the amplification corresponding to P\ =60 deg 
(oblique shock) agrees closer with the measurements than that corresponding to the normal 
shock. 

Numerical simulations by Manning and Lele [26] for a single shear layer-compression wave 
interaction suggest that the radiated acoustic pressure amplitude approaches a saturation level at 
increased distances from the origin of the shear layer (Fig. 7c). The trends of the present model 
qualitatively agree with these results, although in the exponential region, the amplification rates 
appear to be comparatively high with respect to those in the present model. However, a direct 
quantitative comparison of this result with the present model appears difficult in view of the fact 
that a single compression wave is considered in the computation. 

B. Screech Amplitude 

Fig. 8 presents the variation of the excess screech amplitude as a function of Mj as provided 

by the foregoing theory. Experimental data reported by various investigators [10, 40, 43-46] are 
also displayed here in an effort to validate the physical model. The theory for the normal shock 
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case (/?i = 90 deg) suggests that the excess screech amplitude increases from about 17 dB at Mj 
= 1.05 to about 46 dB at Mj = 2.0, with the slope of the curve slowly diminishing with 
increasing A/ y - . The predictions for/?j =60 deg, addressing the sensitivity of the results to the 

shock wave angle, show about 7 dB reduction in the screech amplitude relative to the normal 
shock case. The theory seems to be in reasonable agreement with the data from various sources 
over a range of Mach number, with the normal shock providing the upper limit, and the oblique 
shock (/?i = 60 deg ) yielding the lower limit. 

Measurements for absolute screech amplitude (which includes broadband noise) for a 5:1 
aspect ratio rectangular nozzle obtained by Panda et al. [19] for a wide span of M y confirm the 

general character of the present theory (Fig. 9) in the sense that the screech amplitude increases 
with the jet Mach number until saturation. However, this report does not present the excess tonal 
noise for a direct (quantitative) comparison with the present calculations which provide only the 
excess screech magnitude. It is believed that the saturation (leveling) of the amplitude at high 
Mach number (say M > 1.4), as indicated by the data, may be qualitatively explained to some 
extent if the shock-wave angle varies (decreases) with downstream distance. The relatively rapid 
decline in the amplitude at Mj = 1.7 is perhaps related to the change in the screech mode. 

It is generally believed that the screech tone noises are generated by the interaction between 
the growing instability waves (vortices) developed from the nozzle lip area and the shock cell 
structures. Recently large eddy simulation of jet screech by Berland et al. [29] provided evidence 
to the shock leakage theory proposed in Suzuki and Lele [47]. In the present model, only the 
acoustic waves (with finite amplitude), which is a minor part of the strong hydrodynamic process 
of vortex/shock cell interaction, are extracted, and the shock cells (shocks and expansions) are 
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simplified to normal/oblique shocks. Thus the physics described in the model is quite different 
from the physics of shock/instability wave interaction. In both the situations, the acoustic wave 
intensities somehow grow in an exponential manner before the saturation state. The ability of the 
present model to reasonably predict screech amplitude is perhaps connected with this 
circumstance. In the proposed model, it is thus important to adjust the constants (such as 
M r = 1 AM j in eq. (2b), and the angle (i \ ) to match the experimental data. 

In reality the strongest interaction usually occurs at the 3 rd -5 ,h shock cell, and the shock cells 
become highly deformed and are unsteady. The assumption in the model that these shock cells 
still have the same shape and size and steady is somewhat untenable. 

Also as the acoustic waves are of finite amplitude, the linear theory may not strictly apply. 
Generally speaking, the accuracy of the linear theory appears questionable for weak shocks, 
strong disturbances, and incident angles close to the critical angle [39]. With regard to the 
incident waves of finite amplitude, numerical Euler simulations of Zong et al. [39] with 
amplitudes ranging from 0.01 to 0.5 suggest that the linear theory surprisingly remains reliable to 
extraordinarily large amplitudes. 

The effect of screech mode on the screech intensity is beyond the scope of the present 
investigation. The dependence of screech intensity on jet temperature is also excluded from 
consideration in the present report. A detailed account of nozzle lip thickness is also not 
addressed here. 
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C. Screech Directivity 

Sample calculations of the directivity of screech for a typical value of M j =1.4 for a three- 

source system are depicted in Figs. 10a and 10b, respectively, for the fundamental and the 
second harmonic. An effective value of M c = 0.7 is chosen as a representative value. Results for 

variable strength as obtained from Eq. (10) with a 2 =2.75,03 =5.79 are compared with the 
case of uniform source strength (a 2 = a 3 =1). The values of a 2 and o 3 are obtained from Fig. 6 
at Mj =1.4, relative to the first shock-cell strength of 3.5. For comparison purposes the 
directivity data of Norum [12] for a choked nozzle at Mj = 1.49 are included. The frequencies of 

the fundamental and the second harmonic in this data are respectively 19 Hz and 38 Hz. 

For the fundamental (Fig. 10a), the directivity pattern from the model displays three lobes. 
One lobe, namely the principal lobe, radiates upstream and forms part of the feedback loop. The 
next dominant lobe radiates downstream at a relatively small angle from the downstream flow 
direction. The third lobe is the weakest one, which radiates in a direction normal to the jet axis. 
The directivity of the fundamental displaying the predominant lobe in the upstream direction is 
consistent with the well-known trends of the screech data [12], The relatively rapid decline in 
SPL data at about 20 deg and 160 deg is believed to be a consequence of possible refraction 
effects, which are important for high temperature and high speed jets but are outside the scope of 
the present theory. The data for the fundamental however do not display the weakest lobe 
(normal to the jet axis), as indicated by the theory. It is evident from the present model that the 
effect of variable strength relative to uniform strength is displayed only in the weakest lobe and 
is not significant in the primary lobe directed upstream. 
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Fig. 10b for the directivity of the second harmonic suggests that the peak radiation due to 
screech is in a direction normal to the jet axis — an established result that is in concurrence with 
the available measurements, such as those of Norum [12]. Again, the variation of strength has 
some influence on the predicted directionality of the second harmonic in the three subordinate 
lobes only and does not alter appreciably the directivity of the primary lobe directed normal to 
the jet axis. The theory considerably underpredicts the upstream lobe levels for the harmonic. As 
in the fundamental, the directivity data for the harmonic show possible refraction effects at about 
20 deg and 160 deg. Larger number of sources may in general improve the accuracy of the 
predictions for the directivity- both the fundamental and the second harmonic. 

The proposed method is applicable also to a forward flight stream surrounding the jet by 
appropriately modifying Eq. (2a) and (2b), and with a proper choice for the value of M c in Eq. 

(10). Measurements by Krothapalli et al. [22] with a freestream Mach number from 0 to 0.32 
with an underexpanded heated supersonic jet (nozzle pressure ratios up to 4.5) reveal that the 
screech intensity is unaltered by the forward flight. Cold jet data by Norum and Shearin [48] 
indicate the existence of intense screech tones in simulated forward flight up to Mach 0.41. 

V. Conclusions 

Quantitative estimates for the excess screech tone amplitude over a range of fully expanded 
jet Mach number, as derived from the present model, seem to satisfactorily agree with the 
experimental data for underexpanded supersonic jets. The predicted streamwise amplification of 
the acoustic source strength is also found to be consistent with the available measurements, 
especially at low Mach numbers. This favorable comparison between the theory and the 
measurements suggests that the proposed physical model for the amplification of downstream- 
propagating disturbances of finite amplitude based on shock-acoustic wave interaction theory 
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represents a plausible mechanism of screech in the feedback cycle. The results also indicate that 
the directivity of screech due to variable source strength of the discrete monopole source in each 
shock-cell does not appreciably differ from that based on uniform source strength insofar as the 
principal lobes are concerned. The proposed hypothesis for the determination of screech 
amplitude appears to hold promise in our effort to refine the jet screech modeling. 
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Captions to Figures 


Fig. 1 A typical narrowband farfield shock noise spectrum (adapted from Seiner [1 1]. 

Fig. 2 Schematic of the feedback loop in supersonic jet screech according to Powell [1-3]. 

Fig. 3 Idealized distribution of acoustic sources in an underexpanded jet. 

Fig. 4 Interaction of an acoustic wave with a stationary shock wave (McKenzie & Westphal [34]. 
Fig. 5 Acoustic pressure amplification as a function of upstream Mach number. 

Fig. 6a Streamwise amplification of acoustic source strength from present theory for P\ = 90 deg. 
Fig. 6b Streamwise amplification of acoustic source strength from present theory for (3\ = 60 deg. 

Fig. 6c Comparison of predicted streamwise amplification of acoustic source strength with test 
data. 

Fig. 7a Streamwise amplification of acoustic source strength from the data of Kaji & Nishijima 
[10] for Mj =1.64. 

Fig. 7b Streamwise amplification of acoustic source strength from the data of Panda [19]. 

Fig. 7c Streamwise amplification of acoustic source strength from the numerical simulation of 
Manning & Lele [26], 

Fig. 8 Comparison of the dependence of excess screech amplitude on Mach number. 

Fig. 9 Dependence of measured screech amplitude on Mach number, from data of Panda [19]. 

Fig. 10a Directivity of the screech fundamental. 

Fig. 10b Directivity of screech second harmonic. 
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Fig. 1 A typical narrowband farfield shock noise spectrum (adapted from 
Seiner [11]. 


nozzle 



Fig. 2 Schematic of the feedback loop in supersonic jet screech according to Powell [1-3]. 



Fig. 3 Idealized distribution of acoustic sources in an underexpanded jet. 


refracted sound wave 



Fig. 4 Interaction of an acoustic wave with a stationary shock wave (McKenzie and Westphal 
[34]. 



Fig. 5 Acoustic pressure amplification as a function of upstream Mach number. 
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Fig. 6a Streamwise amplification of acoustic-source strength from present theory for 
P\ = 90 deg. 
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Fig. 6b Streamwise amplification of acoustic-source strength from present theory for 
P\ =60 deg. 


30 



Fig. 6c Comparison of predicted streamwise amplification of acoustic source strength with test 
data. 



FIG. 7a Streamwise amplification of acoustic source strength from the data of Kaji and 
Nishijima for M t = 1 . 64 . 
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Fig. 7b Streamwise amplification of acoustic-source strength from the data of Panda [19]. 
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Fig. 7c Streamwise amplification of acoustic-source strength from the numerical simulation of 
Manning and Lele [26]. 
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Fig. 8. Comparison of the dependence of excess screech amplitude on Mach 
Number. 



Fig. 9 Dependence of measured screech amplitude on Mach number, from data of Panda [19]. 
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Fig. 10a Directivity of the screech fundamental. 
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Fig. 10b Directivity of the screech second harmonic. 
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